A viability criterion for modified gravity with an extra force 
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A recently proposed theory of modified gravity with an exphcit "anomalous" coupling of the 
Ricci curvature to matter is discussed, and an inequality is derived which expresses a necessary and 
sufficient condition to avoid the notorius Dolgov-Kawasaki instability. 
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INTRODUCTION 



modifications of gravity at cosmological 



Recently, moaitications ot gravity at C( 
scales have received much attention jl|, 0, 0] in order to 
explain the cosmic acceleration discovered in 1998 using 
type la supernovae [4]. The alternative is to resort to a 
mysterious form of dark energy with exotic properties ^ : 
a negative pressure P and an energy density p satisfying 
P sa —p and perhaps even P < —p (phantom energy [^). 
The latter easily leads to a Big Rip singularity at a fi- 
nite future [7] . Modified gravity allows one to avoid such 
exotica. The Einstein-Hilbert action 1361 
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is modified to [l|, Q 
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m 

2k 



where /(i?) is an (a priori) arbitrary function of R, 
and the modifications are designed to affect cosmolog- 
ical scales and stay tiny at smaller scales in order not 
to violate the Solar System constraints [8]. The proto- 
type f{R) = R- p*/R (with p Hq 10-33 eV) is 
now regarded as an unviable toy model at best because 
it is subject to a violent instability Q and violates the 
experimental constraints (lol |. 

In order to be viable, modified gravity theories must 
be free of short time scale instabilities and ghosts [ll|, 
12, have a well-posed Cauchy problem 14 1, and have 
the correct cosmological dynamics including an early in- 
flationary era followed by a radiation era, a matter era, 
and a late accelerated era (in many models there are 



problems with the exit from the radiation era 

Modified f{R) gravity comes in three versions: metric 
formalism, in which the action ([2]) is varied with respect 
to the (inverse) metric tensor g""^; Palatini f{R) grav- 
ity, in which variation is with respect to both g"''' and an 
independent, non-metric, connection F^^ but the mat- 
ter part of the action does not depend on F^^ [l^; and 
metric-affine gravity, in which also Cm depends on the 
non-metric connection ^17| . Here we focus on the metric 
approach. 



Recently, Bertolami, Bohmer, Harko, and Lobo (here- 
after BBHL) [111 put a new twist on /(i?) gravity by 
considering the action 
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where fi^2{R) are arbitrary functions of the Ricci curva- 
ture and A is a small parameter (from now on we follow 
(T^ l and set k = SttG = 1). The novelty consists of the 
coupling function f2{R) which adds extra freedom and 
new features. The field equations are 
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where a prime denotes differentiation with respect to R, 

-a 



(2) □ ^ .g-'^V.V,, and T^' = ^ 



Because of 



the extra explicit coupling to matter A/2(i?), T'^™'* is not 
covariantly conserved and energy is exchanged between 
ordinary matter (T^™'') and the "effective matter" repre- 



sented by terms in MR) ™ ^q. (|4|). T^™'' obeys [l 
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The BBHL theory contains intriguing phenomenology: 
all massive particles are subject to an extra force, similar 
to the one arising in scalar-tensor (ST) gravity following 
a conformal transformation to the Einstein frame 1^, 2^ ■ 
In Einstein frame ST gravity the extra force is due to an 
"anomalous" coupling of the matter Lagrangian to the 
Brans-Dicke-like scalar 6, 
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There is, however, an important difference between Ein- 
stein frame ST gravity and BBHL theory: while in the 
former the units of time, length, and mass are not con- 
stant but scale with appropriate powers of the confor- 
mal factor of the conformal trarisformation defining the 
Einstein frame (as explained in [l9| and discussed exten- 
sively in jlH), in the latter there is no such scaling of 



2 



units. For this reason, the BBHL theory [l8| can not 
be reduced to "ordinary" ST or string gravity (in this 
respect, string theory has the same phenomenology of 
ST gravity — indeed, the low-energy limit of the bosonic 
string is a Brans-Dicke theory with Brans-Dicke param- 
eter oj = —1 



where 



2^) 



The extra force on massive particles generated by the 
^f2iR) coupling is always present and causes a devia- 
tion from geodesic paths; therefore, massive test parti- 
cles simply do not exist. Due to this extra force, the 
acceleration law in the weak-field limit of BBHL theory 
assumes a form similar to the one of Modified Newtonian 
Dynamics (MOND) [l^l, which was originally proposed 
to explain galactic rotation curves without dark matter. 
MOND has recently received a relativistic formulation 
in the rather complicated Tensor- Vector-Scalar (TeVeS) 
theory of (H- The BBHL proposal exhibits MOND-like 
phenomenology but has a simpler formal structure than 
TeVeS: as shown below, it amounts to a ST theory with 
two coupling functions, one of which is the coupling of the 
scalar degree of freedom = i? to matter (unorthodox 
in ST gravity [25|). Of course, in order to be viable, the 
BBHL theory must pass the tests mentioned above for 
f{R) gravity and it is not clear yet whether this is possi- 
ble. In this paper we study one of these criteria, namely 
the stability of the theory with respect to local perturba- 
tions. In pure f{R) gravity, a fatal instability develops 
on time scales ~ IQ-^*^ s when /"(i?) < 0. This 
instability, which we refer to as the "Dolgov-Kawasaki 
phenomenon" , was discovered in the prototype model 
f{R) = R — ^i'^/R which is ruled out (and only cured 
by adding extra terms to f{R) [U, 12, (111), and then 
generalized to arbitrary f{R) models 27|. For the BBHL 
theory, the corresponding stability criterion turns out to 
be f'^(K) + 2Xf^{R) > (see Sec. 3; see Refs. [28] for 
other types of instabilities). 



EQUIVALENCE WITH AN ANOMALOUS ST 
THEORY 



It is well known that pure f{R) gravity ([T]) is equivalent 
to a ST theory [2^; here we revisit this equivalence and 
generalize it to BBHL theory. 

By introducing a new field = (j), the action ([3]) is writ- 
ten as 

■0) + [i + A/2(0)]/:„ 

(7)' 

and, further introducing the field "0(0) = fi{4') (where 
now a prime denotes differentiation with respect to (f> 
[s^]), one can write 
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with (/)(■(/') given by inverting ipi4') = /((</')■ The ac- 
tions ([3|) and ([8|) are equivalent when /"(i?) ^ 0: in fact, 
by setting (j) = R, eq. ^ reduces trivially to eq. ([3]). 
Vice-versa, variation of ([7]) with respect to (j> yields 



(i? -(/.)/{' (0) + 2A/^((^)/:,„ = O 
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In vacuo (Cm = 0), this equation yields (j) — R when- 
ever /" ^ 29] . In the presence of matter there seem 
to be other possibilities which are, however, excluded as 
follows. When £,„ ^ 0, the actions ([3]) and ([T]) are equiv- 
alent if [R - (f)) f['{(j)) + 2 A/^' ((/))£„ ^ 0. When eq. ^ 
is satisfied, wc have a pathological case which, upon in- 
tegration of this equation, corresponds to 
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But if eq. holds, then the action ([7]) reduces to pure 
matter without the gravity sector and we dismiss this 
case. Then, the actions ([3]) and ([U are equivalent when 
f['(R) ^ 0, as in pure /(i?) gravity [SiLThe action ^ 
corresponds to a Brans-Dicke theory [30| with a single 
scalar field, vanishing Brans-Dicke parameter w, and an 
unorthodox coupling U{tjj) to matter. Actions of this 
kind have been contemplated before 31 , 13, 2^ , but little 
is known about them. 



BBHL THEORY AND INSTABILITIES 

The trace of the field equations is, in terms of i?, 

3 [f['{R) + 2XCraf'^{R)] UR + 3 [f{\R) + 2AA„/^'(i?)] 
V'RVcR + 12X f^{R)V^C„y,,R + f[{R)R - 2/i(i?) 
+2A£„/^(i?)i? = [1 -I- Xf2{R)] T(") ~ 6A/^(i?)n£„ , 

(13) 

where T^™' = T('")''a ■ As customary in f{R) gravity, 
we parametrize the function fi{R) as fi{R) = R+eip{R), 
where e and A must necessarily be small to respect the 
Solar System constraints [11]. Following 0, we expand 
the spacetime quantities of interest as the sum of a back- 
ground with constant curvature and a small perturba- 
tion: R = Rq + Ri, T ~ To + Ti, Cm = Co + Ci, 
and the spacetime metric can locally be approximated by 
dab — fjab + habj whcrc Tjab is the Minkowski metric. There 
are really two approximations here. The first is an adi- 
abatic expansion around a de Sitter space with constant 
curvature, which is justified on timescales much shorter 
than the Hubble time. The second is a local expansion 
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over small spacetime regions that are locally flat (hence 
the appearance of rjab)- These approximations are com- 
mon in f{R) gravity {e.g., [1,[23) ^^id 1980s literature 
on inflation. Accordingly, /i(i?) = Rq + Ri + eip{Ro) + 
t^'{Ro)Ri + f[{K) = 1 + e(^'(i?o) + eip"iRo)Ri + ... 
and the linearized version of the trace equation in 
the perturbations becomes 

3 [e^"iRo) + 2A/^'(i?fl)] DR^ + [eip"{Ro)Ro - 1 
-e^'iRo) + 2Xf!,{Ro)Co + 2XCof^{Ro)Ro 
-Xf^{Ro)To + 6A/^'(Eo) {aCo)] Ri 
= -2Xf^{Ro)RoCi + [1 + XMRo)] Ti 
-6Xf!,iRo)DCi , (14) 

where the zero order equation 

f[{Ro)Ro - 2/i(i?o) + 2XCof^iRo)Ro = [1 + A/2(i?o)] Tq 

(15) 

has been used. Eq. (fTi|) is further rewritten as 

Ri - V^Ri + mlffRi = {3 [e(^"(i?o) + 2A/^'(i?o)]}"' 
{2Xf!,[Ro)Ci - [1 + A/2(i?o)] Ti + 6A/^(i?o)n/:i} (16) 

where the effective mass We// of the dynamical degree 
of freedom Ri is given by 

miff = {3 [e^"(i?o) + 2A/^'(i?o)]}'' {1 + e(^'(i?o) 
+e(^"(i?o)i?o - 2XCa [f^^iRo) + (^o)i?o] + A/^(i?o)ro} . 

The dominant term on the right hand side is 
{3[e(^"(i?o) + 2A/^'(i?o)]}"^ and the effective mass 
squared must be non-negative for stability. Therefore, 
eip"{R) + 2A/2 (i?) > is the stability criterion for the 
BBHL theory against Dolgov-Kawasaki instabilities. 

OUTLOOKS 

The inequality eip"{R) + 2A/2'(i?) > generalizes the 
stability condition f "(R) = eip"{R) > found in pure 
f{R) gravity The survival of BBHL theory 

is subject to satisfying the other (physically independent) 
viability criteria mentioned above, which require a sepa- 
rate analysis and will be analyzed in future publications. 
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